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We derive the high-temperature expansion of the Helmholtz free energy up to the order /3 17 of the one- 
dimensional spin-S Ising model, with single-ion anisotropy term, in the presence of a longitudinal magnetic 
field. We show that the values of some thermodynamical functions for the ferromagnetic models, in the pres- 
ence of a weak magnetic field, are not small corrections to their values with h = 0. This model with S = 3 
was applied by Kishine et al. [J.-i. Kishine et al., Phys. Rev. B, 2006, 74, 224419] to analyze experimental data 
of the single-chain magnet [Mn(saItmen)] 2 [Ni(pac) 2 (py)2](PF6)2 for T < 40 K. We show that for T < 35 K 
the thermodynamic functions of the large-spin limit model are poor approximations to their analogous spin-3 
functions. 

Key words: quantum statistical mechanics, one-dimensional Ising model, spin-S models, large-spin limit Ising 
model, single-chain magnets. 
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1. Introduction 

The one-dimensional spin-1/2 Ising model with first-neighbor interaction, in the presence of a lon- 
gitudinal magnetic field, was exactly solved in 1925 |1]. The Helmholtz free energy (HFE) of this model 
with S — 1/2 has a simple mathematical expression |2|. The exact expression of the HFE for the S — 1 
ferromagnetic model in the presence of a longitudinal magnetic field was derived in 1976 by Krinsky and 
Furman [ 3] using the matrix density approach. (This work seems to have been neglected by the subse- 
quent literature, though). More recently, the HFE of the S = 1 flU] and S = 3/2 @] of the Ising model in 
the presence of an external magnetic field have been written as a set of coupled equations and solved 
numerically. 

In 2007 Rojas et al. [7] published the high-temperature expansion of the HFE of the Ising model for 
arbitrary value S of the spin, in the absence of a magnetic field, up to order /3 40 . The nice feature of such 
expansion is that the value S of the spin is an arbitrary parameter (S = 1/2, 1, 3/2,...). The absence of 
an external magnetic field in the model discussed in [7] yields no distinction in the behavior of some 
thermodynamical functions between the ferromagnetic and the anti-ferromagnetic (AF) case, e.g. the 
specific heat per site. 

An active area of molecular chemistry is that of designing new magnetic materials that present a one- 
dimensional nanomagnetic behaviour with a strong anisotropy axis (§]. Some molecular single-chain 
magnets (SCMs) exhibit a strong uniaxial (Ising) anisotropy. Kishine et al. @] applied the Blume-Capel 
model (iji EH (the one-dimensional Ising model with the single-ion anisotropy term) with S = 3 in 
the presence of a longitudinal magnetic field to analyze the low-energy dynamics response of the SCM 
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[Mn(saltmen)]2[Ni(pac)2(py)2](PF6)2 for temperature T < 40 K. They concluded that the experimental 
data support the view that this SCM can be described, in this window of temperature, by the spin-3 of this 
model. 

We apply the method of cumulants described in [12] to calculate the high-temperature expansion of 
the HFE of any one-dimensional Hamiltonian which is invariant under space translation and satisfies 
the periodic space condition. In this approach, in order to obtain the exact coefficient that multiplies the 
term of order p" in the expansion we have to calculate a set of functions named HrL, m- 1,2,..., n. The 
interested reader will find a survey of the method in reference 11311 . In section [2] we apply the results of 
reference (l^] to calculate the high-temperature expansion of the HFE of the normalized one-dimensional 
spin-S Ising model with single-ion anisotropy term in the presence of an external longitudinal magnetic 
field up to order B 17 . The nice feature about this expansion is that it is faster to apply than the transfer- 
matrix method |2], although the latter one provides exact curves in the whole interval of temperatures. 
In order to show the importance of calculating the high-temperature expansion of the spin-S Ising model 
in the presence of an external magnetic field, in section [3] we compare the behavior of certain thermo- 
dynamical functions of the ferromagnetic and of the AF Blume-Capel models in the presence of a weak 
longitudinal magnetic field. This comparison is made for various values of the spin, including the large- 
spin limit (S — >■ oo). Section|4]presents the thermodynamics of the SCM (for spin-3 and the large-spin limit 
model) in the interval 11 K < T < 40 K. Finally, in section[5]we present our conclusions. 



2. Thermodynamics of the spin-S Ising model with a single-ion 
anisotropy term 

The Hamiltonian of the spin-S Ising model with a single-ion anisotropy term in the presence of a 
longitudinal magnetic field is |13]] 

N 

H' S =E 7'SfS? +1 -fe'Sf + D'(Sf) 2 , (1) 

;' = 1 

where S? is the z component of the spin S operator with norm: ||S|| 2 = S(S+ 1), S - \, 1, |,2, at the 
i-fh site of the chain; /' is the exchange strength and it can have negative value (ferromagnetic model) or 
positive value (AF model). In references Hill and Hill we studied the exact thermodynamics of this model 
for S = | and 1 (with h = 0), respectively. 

If the large-spin limit (S — ► oo) is applied directly to the Hamiltonian Q}, all the thermodynamic func- 
tions of the large-spin limit model will diverge. In order to keep the functions finite in this limit, we study 
the normalized version of the Hamiltonian CD, that is, 

jv 13 
Hs=E \j s i s li- hs i +D ( s if}' S=-,l,-,2,..., (2) 

where s? is the z component of the spin operator ? that has norm 1. The s operator is denned as 

S 1 3 

S=-,l,-,2,.... (3) 



vW+T) 2 2 

Making S -«■ oo in the Hamiltonian CD we obtain its large-spin limit. 
The Hamiltonians CD and {2} are identical when 



/=S(S+1)/', h = v'SiS+Dti and D = S(S+1)D' (4) 



withS= i,l,§,2,.... 



Let Z' s (/', h', D'; p) be the partition function derived from Hamiltonian CD, Z' s (/', h', D'\ 5) = Tr (e _/3H 's j . 
Its HFE, in the thermodynamic limit, is called W' s , where 

W^f,ti,D';P) = - lim ^ \ hx[Z' s {j',h',D';P)\. (5) 

N^oo iV p 
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The analogous functions for the normalized Hamiltonian l|2} are 2 S (J, h,D; p) = Tr(e ^ Hl ), and 

W s {J,h,D\p) = - lim ^ \ ln[2 s U,h,D;P)] . (6) 

iV— oo JM p 

For arbitrary spin value S, the relation between the HFE's (5} and GD is 

%U',h',D';P) = W s (S{S+l)j', \/S{S+l)h',S{S+l)D';P). (7) 

The Hamiltonian (1) and its normalized version l|2) belong to a subclass of chain Hamiltonians that 
satisfy periodic conditions and that can be decomposed as H^+i = P, ii+ i + Q; , in which the operator 
Pj,i+i depends on two sites (the i-th and the (i + l)-th sites), the operator Q ; - depends only on the i-th site; 
moreover, 



tr,(P w+ i)=0. (8) 



It is simple to show that the functions H[ n f } n calculated in the method presented in references Il2l.ll; 
for the subclass of Hamiltonians satisfying the condition l(8) can be written as 

rj(n) _ rr(l) rr(n-l) , V- / n 12 n 23 m,m+l \ . . 

H i, m - H i,i xH i,m-i + L \7rr~^j ••• w i / (9) 

ni+W2+rt3+...+n m =n 
n,#0, i=l,..., m 



with 2« m«H, 

Equation (9) tells us that for Hamiltonians satisfying condition (§)> the function H^ m has a contribu- 
tion from ff^H , which are calculated to obtain the coefficient of order in the expansion of the HFE. 
The number of terms to be calculated in each order of /3 in this expansion is then remarkably reduced. 

The result (9} permits us to write: 

H["l=(H$) n , n= 1,2,3... (10) 

and 

H i"Li = H iA x<„"- 1 2 ! + ^(H? 2 H23...H„_ 1 ,„) g> (11) 
where the g-trace ((. . .} g ) in equations (9) and dTT) means 

"—'))' (12) 

where £™j = « and + with j = 1,2,..., ra. The notation 3 a (H" 1 2 ,H" 2 3 ,...,H^" m+1 ) means all the 
distinct permutations of the n operators where m of them, {Hi 2,H2,3, . . . ,H m>m+ i}, are different. The 
notation (...) corresponds to calculating the normalized traces on the indexes: \,2,...,m+ \ 

The results l[9)-l[ll} are valid for any ID Hamiltonian that is invariant under space translation, satisfies 
the periodic space condition and the condition (8). 

In this article we study a ID model that satisfies the necessary conditions to apply the method of 
reference 1 12] plus the condition 1(8). All the terms in this Hamiltonian are commutative; hence, the g- 
traces in H^-fTO can be replaced by the usual normalized traces 11211 . 

We calculate the high-temperature expansion of W s {J,h,D;p), for arbitrary value of the spin S , up 
to order /3 17 . The relation l[7) can be applied to yield the high-temperature expansion of W^U' , h',D'; p) 
from the expansion of W S {J, h, D; /3). 
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Here we present the high-temperature expansion ofW s {J, h,D;/3), for any spin value S,up to order /3 2 , 



ln(2S+l) D 
W s (J,h,D;P) = ( — - + - + 

t2 i n2 



1 Sh 2 1 D 2 2 SD 2 
■ + ■ 



6 S+l 30 S(S+1) 45 S+l 



i r id 2 i s/ 2 i 



18 S+l 45 S+l 18 S+l 6 S+l 



1 J 2 D 



135 (S+l) 2 



1 h 2 D 2 S 2 /i 2 L> 4 Sfc 2 D 2 /Sfe 2 1 ] 2 D 

■ + T + T + - 



30S(S+ 1) 2 45 (S+l) 2 45 (S+l) 2 9 (S+l) 2 45 S(S+1) 2 
1 h 2 D 8 SJ 2 D 4 D 3 1 /fc 2 4 S 2 / 2 D 



■ + - 



90 (S+l) 2 135 (S+l) 2 405 (S+l) 2 9 (S+l) 2 135 (S+l) 2 

1 D 3 1 JS 2 h 2 4 D 3 16 SD 3 

126 S 2 (S+1) 2 9 (S+l) 2 315 S(S+ 1) 2 2835 (S+l) 2 
S 2 -D 3 1 

P 2 + &(P 3 ). (13) 



2835 (S+l) 2 

We should note that this high-temperature expansion is valid for positive, null or negative values of 
/, and for arbitrary values of S, h and D. The HFE of the large-spin limit model of Hamiltonian {2} is 
calculated from the high-temperature expansion of W s (/, h, D; p) by taking the limit S — «■ oo. 

The authors maintain a websit^l in which the interested reader may find data files on the arbitrary 
finite spin-S and the large-spin limit (S -► oo) HFE's of the normalized Hamiltonian (2} up to order /3 17 . 

We have a few general comments on the function W S {J, h, D; /3): 
the expansion (13\ of the HFE is an even function of the external longitudinal magnetic field h; 
ii) the function W s {], h, D; (i) is even in the parameter / for an external magnetic field with null longitudi- 
nal component (h = 0); 

Hi) even for h-0, the HFE (13) is sensitive to the sign of the parameter D. 

By direct comparison, we verify that our expansion of the HFE of the spin-S Ising model coincides for 
S = 1/2 and S = 1 with the high-temperature expansions of the exact results of references (3 and (H, 
respectively. 

It is simple to understand why the second comment is valid for the exact expression of the HFE of the 
model {2). The partition function comes from the calculation of the traces of operators (H s )". In refer- 
ence llill we showed that for any spin-S only tr,- (s?) 2 ' + 0. In the absence of the longitudinal component 
of the magnetic field {h - 0), only products with even number of operator /-s?s? +1 give non-null contribu- 
tions tothetr[(H s )"]. 

The contribution of the single-ion anisotropy term in Hamiltonian {2) to the partition function comes 
from the operator exp{ - [/3D[S^) 2 ]/[S(S+ 1)]}. For positive values of the crystal field (D > 0) the main 
contribution of this operator to the partition function comes from the smallest eigenvalues of (S?) 2 . On 
the other hand, for negative values of D the states with the largest eigenvalues of (S?) 2 are favored. Such 
distinct behavior for positive and negative values of D explains the origin of condition (Hi). 



3. The ferromagnetic and anti-ferromagnetic models in the presence of 
a weak magnetic field 

The exact expression of the HFE of the Hamiltonian (2} and its thermodynamic functions are un- 
known for an arbitrary value of spin 5. Our work has been that of calculating the high-temperature 
expansion of thermodynamic functions of one-dimensional models. In section [2] we mentioned that the 
HFE's of the Hamiltonians (1} and (2) in a magnetic field with null longitudinal component {h - 0), is an 
even function of /. As a consequence, for h = 0, several thermodynamic functions are the same for the 
ferromagnetic (/ < 0) and AF (/ > 0) models, namely: the specific heat, the internal energy, the entropy 
and the mean value of the square of the z component of the spin ((S z ) 2 ). For h = 0, we also have that the 

1 http://www.proac.uff.br/mtt 
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correlation function between first neighbors satisfies the equality 

<€ s {.-J,Q,D;p) = -<€ s U,0,D;P) (14) 

for s= |,1,§>2,.... 

The high-temperature expansion (13) of the HFE of the spin-s model is valid for arbitrary values of h. 
In order to verify how important the presence of a longitudinal magnetic field is to the thermodynamic 
properties of the ferro and AF models, in this section we examine how the thermodynamic functions 
of the ferromagnetic and AF models differ from their respective values at h = when they are in the 
presence of a weak magnetic field for different values of spin. 

Throughout this section we consider / = 1 (AF model) or / = -1 (ferromagnetic model). Again, the 
parameters are in units of |/| and the expansions are in powers of 

Let & s {J,h,D; /3) be a given thermodynamic function with spin-s (see relation {3)) derived from 
Hamiltonian (2). Its percentage difference to the value of & S {\,Q,D; j6) is defined as 



& s (l,0, D;p)-& s U,h,D;p) 



Let us compare the thermodynamic functions. 



x 100%, /=+l. (IS) 



1) Comparison of the specific heat per site: C S [J, h,D; f3). 

For each spin-s, the specific heat functions of the ferromagnetic and AF models, in the high temper- 
ature limit (/3 — ► 0), reach the same value for any h and D. The p 2 term in this thermodynamic function 
has a J 2 dependence for both models. 

Figure[l]show AC s (+l,ft/|/|,D/|/|; |/|/f) for s = 1/2,2,4 and oo (large-spin limit model) as a function 
of hl\J\ and (|/|jS). In figure Q](a) we have D = 0, \J\/3 = 1.6 and hl\J\ e [0,0.1]. For these spin values, 
|AC S (-1, ft/1/1,0; 1.6) <43% for the ferromagnetic models, whereas |AQ(l,ft/|/|,0;1.6)| < 0.42% for the 
AF models. In the AF case the function AC S (1, h/\J\,0; 1.6) for s = 1/2 very closely approximates its large- 
spin limit (s — >■ oo) version. Figure[l](b) shows the percentage difference AC 5 as a function of (|/|/3), for 
DI\J\ = -0.5 and hl\J\ - 0.1. Again we have that the percentage differences of the ferromagnetic models 




\J\ \J\ P 

Figure 1. The percentage difference of the specific heat for h = and for a weak magnetic field h for the 
AF (/ = 1) and for the ferromagnetic models (/ = -1). Notice that for the AF case, all curves are close 
to the horizontal axis. This comparison is made for both models with s = 1/2 (dashed line), 5 = 2 (dash- 
dotted line), 5 = 4 (solid line) and the large-spin limit model (dotted line). In figure[l](a) this percentage 
difference is plotted as a function of hl\J\ with D = and (|/|/3) = 1.6. Figure[l](b) presents the curves of 
AC S versus (|/|j6) with DI\J\ = -0.5 and hl\J\ = 0.1. 
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are high (| AC S (- 1,0.1, -0.5; <38% whereas the percentage differences for the AF models are much 
smaller (|AC S (1,0.1, -0.5; |/|j8)| <3%). 

The previous discussion exemplifies the fact that the high-temperature expansion of the specific heat 
of the AF model in the presence of a weak longitudinal ma gnet ic field can be approximated by the corre- 
sponding expression derived from the results of reference 11311 (where we have h - 0). Figure[T]shows us 
that for arbitrary spin-.? we can improve this approximation by recognizing that C s (l, hl\]\,DI\J\; « 
Cs(l,0,DI\J\;\J\P) + ACii 2 il,hl\J\,D/\J\; \J\fi), at least in the region in which hl\J\ « 1 and \J\P<1.6. It 
is important to recall that the exact expression of C1/2U, h, D; /3) is known fl3l . 

For the spin-s ferromagnetic model (2} we need to know the high-temperature expansion of the HFE 
in the presence of longitudinal magnetic field to obtain the value of C s {-\,hl\J\,DI\]\; \J\f5), even in the 
presence of a weak field. 

2) Comparison of the correlation function between first neighbors per site: S^fJ, h,D; 
Equation fl4t tells us that in the absence of magnetic field (h = 0), 

^ s (-/,O,D;/3) = -<if? s (/,O,D;,0) (16) 

for s = |,1,|,2, ... Result fl6t describes the parallel (anti-parallel) alignment of the neighboring spins in 
the ferromagnetic (AF) model. 

We want to check how the function ^(+1, h/\J\,D/\J\;\J\f}) differs from the equality (16) in the pres- 
ence of a weak magnetic field. In order to quantify this deviation, we define 



3}^ s {hl\J\,DI\J\;\J\fi) = 



^ s [-l,hl\J\,D/\J\;\J\p)+^ s [l,h/\J\,DI\J\;\J\p) 
<e s l-l,hl\J\,DI\J\;\J\p) 



x 100%. (17) 



In figure[2](a) we show the percentage function fFJ) . with D = and DI\J\ = -0.5, with h/\J\ = 0.1, to 
verify the departure from equation ( Tl6"l for the ferromagnetic and AF models in the presence of a weak 
magnetic field. The function ( 2i c € s (h I \ J\ , D I \ J\ ; | J\ p) depends on the spin value and for lower temperatures 
(|/|)0 ~ 1.9) it can be around 10%. 

3) Comparison of the magnetic susceptibility per site: XsU, h,D; p). 




< 



(b) 



// / 



0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 2. Left: comparison of the percentage function ©"rfs for the s = 2 and the large-spin limit (s — co) 
models for hl\J\ = 0.1: for s -2, the cases D = (dashed line) and DI\J\ = -0.5 (solid line) are shown, 
whereas for the s — 00 model the D = (dotted line) and DI\J\ - -0.5 (dash-dotted line) cases are shown. 
Right: comparison of the percentage difference of the magnetic susceptibility A% s (-l,0.1,D;/3) in the 
ferromagnetic model (/= -1) for se {2,4,oo}, for DI\J\ = -0.5 (upper set of curves) and for DI\J\ = -0.5 
(lower set of curves); in both sets, the s = 2 (dashed line), .s = 4 (solid line) and s — 00 (dotted line) cases 
are shown. 
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It is simple to obtain the high-temperature expansion of the magnetic susceptibility per site, 
X s {J,h,D;P), from the series (13) for the HFE {% s tp) = -d 2 W s {p)ldh 2 ). 

At h = 0, the magnetic susceptibility of the ferromagnetic and AF models are distinct. We define 
a percentage difference analogous to equation (15) to present in figure [2] (b) the difference between 
j s (-l,0.1,D/|/|;|/|/3) and x s (-l,0,D/|/|; \J\fi) in the ferromagnetic models with 5 = 2,4 and the large- 
spin limit model. In this figure we have DI\J\ - +0.5. For lower temperatures (|/|/3 ~ 2) we have percent- 
age differences around 20%. By looking at figure[2](b) we see that for J - -1 we have d(Aj s )/d/3| D= 5 > 
d(Ax s )/d/3| D=0 5 . This inequality is explained by the fact that for D < the states with the largest values 
of (S? ) 2 are more probable, whereas the states with the smallest values of (S?) 2 are favored in the model 
with D > 0. 

For the AF models (/ = 1), in the interval \J\/3 e [0, 1.43], with hl\J\ = 0.1 and DI\J\ = +0.5, we obtain 
from the percentage difference <fT5> |Aj,fl. 0.1, DI\J\;\J\p)\ <0.6%. 

We point out that to derive the high-temperature expansion of XsU ,h,D;p) from the HFE, we need 
the dependence of the HFE on h. That is not the case of the previous paper on the spin-S Ising model (7||. 



4. Thermodynamic behavior of the single-chain magnet in the interval 
10 K£ T£40 K 

The Hamiltonian (l), with S = 3, was applied by Kishine et al. [9] to analyze the low energy dynamics 
of the single-chain magnet (SCM) [MnfsaltmenJhlNHpao^fpy^KPFeh for T < 40 K in the presence of a 
weak magnetic field. Their Hamiltonian (2) is identical to ours QJ once we set /' = -2t. 

In 11911 we studied the thermodynamics of the unitary spin-s XXZ model with a single-ion anisotropy 
term in the presence of a magnetic field in the z direction. We obtained that the specific heat, the magneti- 
zation and the magnetic susceptibility of this model with s-3 are well approximated by their respective 
large-spin limit versions, in the temperature range of \J\/3 < 1. For the non-normalized S = 3 XXZ model 
this range corresponds to \J\f3 < 0.083 (see section 2.2 of reference 1 19]). 

In this section, we compare the Hamiltonian (I) having spin-3 with its large-spin limit version for T < 
40 K. By "large-spin limit model" we mean that the z component of its spin vector S varies continuously, 
namely, S? = 2 \Z3cost0,), in which 0,- e [0,n], and i e {1,2,..., A/}. In order to relate our analysis to the 
aforementioned SCM we use the same parameter values as in (gj], namely, 

/' 

V = -1.6K (18a) 
k 

and 

D' 

— = -2.5 K, (18b) 
k 

in the spin-3 and large-spin limit Hamiltonians. One is reminded that k is the Boltzmann constant. Taking 
the value (18a) for /', for T ~ 40 K we have \J'\fi ~ 0.04. Note that the temperature region characterized 
by 1/1/3 < 0.04 is contained in the temperature region in which the spin-3 XXZ model behaves very much 
like its large-spin limit version. 

Along this section, we take h' Ik - 0.25 K, which is a weak magnetic field (h'/\D'\ = 0.1). 

Let ^s(/i;/3) be a thermodynamic function. Its percentage difference of the S-3 and the large-spin 
limit models is defined as 



&s~oo(h;p)-& 3 (h;p) 



x 100%. (19) 



In figure [3] we plot the percentage difference (19) of the specific heat, the z component of the mag- 
netization per site M l z S] U' ' ,ti ',£>'; /3) = -^^) and the magnetic susceptibility per site. The 
dash-dotted curve corresponds to the specific heat in the interval T e [14.1 K,40 K]. For T < 40 K or 
1/1/3 > 0.046, we verify that AC(0.25; fi) > 10%. The percentage difference of the z component of the mag- 
netization, AM z (0.25;/3), is given by the solid line in the figure for T e [11.63 K,40 K], while the dotted 
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curve describes Ax (0.25; p) in the temperature window T e [12.99 K,40 K]. We verify that for T < 19.3 K, 
that is, I/' |)3 > 0.083, the z component of the magnetization and the magnetic susceptibility of the spin-3 
model of the SCM differ more than 8.5% from their corresponding large-spin limit values. 




T{K) 

Figure 3. The percentage differences between the thermodynamic functions of the spin-3 model and 
its large-spin limit version of the SCM. The dash-dotted line corresponds to AC(0.25; fi), the solid line to 
M4 Z (0.25; p) and the dotted line to Aj(0.25; fi). 



Due to the ratio \J'\I\D'\ value of 0.64 for this SCM, the single-ion anisotropy term gives the main 
contribution to its thermodynamics. Since the crystal field D' is negative (see equation ( 118b) ). the states 
with S* — ±3 and +2 are the most probable. The modulus of the amount of energy, in units of k, for the 
single-ion anisotropy term to change from states with S? = +2 to states with Sf = +3, and vice-versa, is 
12.5 K. This value is close to one-third of 40 K. We conclude that for T < 40 K (that is, \ J'\p > 0.04), the 
discretized nature of the spin-3 is still an important feature of the model. This result is very different 
from that of reference Hill for the spin-3 XXZ model with a single-ion anisotropy term in the presence 
of a magnetic field in the z direction. 

Our expansion of the HFE with S = 3 can be easily used to fit the experimental data of this SCM. It can 
also be applied to determine the best fit of the parameters /' and D'. 



5. Conclusions 

The method of reference [12,] permits one to calculate the high-temperature expansion of the HFE of 
any chain Hamiltonian with interaction between first neighbors which is invariant under space transla- 
tion and satisfies a periodic space condition. If the Hamiltonian also satisfies the condition 0, we show 
that some of the terms that contribute to the function H?rL have already been calculated in a lower order 

of n, that is, H^Ir?-.. A set of important ID Hamiltonians satisfy the condition 0. This contribution from 
lower order of n to this auxiliary function permits us to compute the HFE of the one-dimension Ising 
model with single-ion anisotropy term in the presence of a longitudinal magnetic field up to order /3 17 , 
for arbitrary values of spin S and of the other parameters /, h and D. Upon performing the numerical 
analysis of the thermodynamics of spin models (1} and one must know the value of the spin before- 
hand. 

We discuss the thermodynamics of the normalized Hamiltonian (see equation 10), but equation 
relates the HFE of the Hamiltonian (I) and its normalized version (equation 0). 

Some thermodynamic functions are insensitive to the sign of / in the absence of a magnetic field; for 
instance, the specific heat of the ferromagnetic (/ < 0) and AF (/ > 0) spin-S Ising models with single-ion 
anisotropy term are identical with h = 0. In the absence of a magnetic field, the ferromagnetic model 
favors parallel neighboring spins while in the AF model the anti-parallel pairs are more probable. The 
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effect of the presence of an external magnetic field in both models is that of favoring the alignment of 
spins at each site to the field direction. In the ferromagnetic model such alignment is favored by the cou- 
pling between neighboring spins. On the other hand, in the AF model there is a competition between 
the coupling of neighboring spins (favoring anti-parallel alignment) and the Zeeman term (forcing all 
the spins in the chain to align with the external magnetic field). As a consequence of this competition, 
there is a perturbative effect in the AF model due to the presence of a weak magnetic field; in such 
regime, their thermodynamic functions ^$(1, hl\J\,hl\J\\\J\p) for the spin-5 model can be approximated 
by & s Q.,hl\J\,DI\J\;\J\p) « & s (\,0,DI\J\;\J\p) + A&u 2 {\,hl\J\,hl\J\;\J\p). The exact expression of the 
HFE of the spin-1/2 Ising model with a single-ion anisotropy term in the presence of a longitudinal mag- 
netic field is known In the ferromagnetic model, the effect of the Zeeman coupling cannot be treated 
anymore as a perturbation to the interaction between first neighbors and to the single-ion anisotropy 
term for >0.04. 

Kishine et al. 190 applied the spin-3 Hamiltonian Q) to analyse the low energy dynamics of the 
[Mn(saltmen)] 2 [Ni(pac) 2 (py)2](PF6)2 SCM for temperatures T < 40 K. Although S - 3 could be consid- 
ered a high spin value in the region of T ~ 40 K (|/'|/3 ~ 0.04), the negative value of D' favors the states 
with S? = +2 and +3. The modulus of the amount of energy required by the single-ion anisotropy term 
to have S z varied from +2 ^ +3 and vice-versa is about one-third of the thermal energy available for 
T < 40 K. Our results show that the large-spin limit model, with the spin-3 replaced by the classical vector 
in the Hamiltonian, yields a poor approximation to the behavior of this SCM for T < 18 K and for the set 
of parameters values <18a> and <18bt . 

Finally, it is very important to point out that our high-temperature expansion of the HFE of the spin-5 
Ising model can be applied to fit experimental data of new materials with one-dimensional behavior and 
strong anisotropy axis. This expansion leaves the spin value of the material as one of the parameters to 
be determined by the best fit. The expansion is valid for positive and negative values of / and D in the 
presence of a longitudinal magnetic field that does not have to be a weak field. 
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Bi/icoKOTeMnepaTypm/iM po3i</iafl fl/ia K/iaci/iMHoY MOfle/ii l3mra 

i3 H/ieHOM S 2 Z 

M.T. ToMaP, 0. PoxacP 

iHCTHTyT c(>i3HKW : ctJeflepa^bHwR yHiBepcwTeT <t>/iyMiHeHce, HiTepofi-RJ, Epa3w/im 
<t>aKy/ibTeT tohhux HayK, ct>eflepa/ibHni7i ymBepcwTeT m. /laBpac, /laBpac-MG, Epa3n/im 

Mn BUBOflWMO BucoKOTeMnepaTypHUM po3K/iafl Bi/ibHo'i eHeprii r e/ibMro/ibL(a flo 4/iei-iiB nopnflKy /3 17 pjw\ op,no- 

BMMipHOI S-CniHOBOl' MOfle/li l3iHTa, i3 OflHOiOHHOK) aHi30TpOni£K) B npi/lCyTHOCTi n03flOB>KHbOrO MarHeTHoro 

no/in. Mn noKa3yeMO, luo 3Ha4eHHfl TepMOflnHaMHHnx 4>yHKi4ii?l (fiepoMarHeTHnx MOfle/iefi b npncymocTi aia- 
6oro MarHeTHoro no/ia He e Ma/inMn nonpaBKaMn npn h = 0. Up MOfle/ib 3 S = 3 6y/ia 3acTOCOBaHa KawiHe 
Ta iH. [J.-i. Kishine et al., Phys. Rev. B, 2006, 74, 224419] fl/w aHa/ii3y ei<cnepnMeHTa/ibHnx flaHnx OflHO/iaH- 
L4KJXKOBOro MarHeia [Mn(saltmen)]2[Ni(pac)2(py)2l(PF6)2 n P w T < 40 K. Mn noKa3yeMO, 1140 npn T < 35 K 
TepM0flHHaMi4Hi cf>yHKi4iV MOfle/ii b rpaHni^i Be/iwKoro cniHy e noram/iM Ha6/in>KeHHflM fl/ia aHa/iori4Hnx flo hux 
3-cniHOBnx cbyH Ki^i 

K/itoHOBi cjioBa: KBamoBa CTaTucTmna Mexaw/ra, o^HOBUMipHa MO^enb l3iHra, S-cniHOBi MOfleni, rpaHULfa 
Beni/iKoro cniHy, OAHOJiaHU,K»KKOBi Mamem 
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